SCIENTIFIC 

REPORTS 





OPEN 



SUBJECT AREAS: 

SURFACES, INTERFACES 
AND THIN FILMS 

STRUCTURAL PROPERTIES 



Received 
19 March 2014 

Accepted 
14 July 2014 

Published 
1 August 2014 



Correspondence and 
requests for materials 
should be addressed to 
C.L. (chenglu@uow. 

edu.au] 



Atomistic Simulation of Tensile 
Deformation Behavior of ^5 Tilt Grain 
Boundaries in Copper Bicrystal 



Liang Zhang, Cheng Lu & Kiet Tieu 



School of Mechanical, Materials and Mechatronic Engineering, University of Wollongong, Wollongong, NSW 2522, Australia. 

Experiments on polycrystalline metallic samples have indicated that Grain boundary (GB) structure can 
affect many material properties related to fracture and plasticity. In this study, atomistic simulations are 
employed to investigate the structures and mechanical behavior of both symmetric and asymmetric 
^5[0 0 1] tut GBs of copper bicrystal. First, the equilibrium GB structures are generated by molecular 
statics simulation at OK. The results show that the ^^5 asymmetric GBs with different inclination angles ((|)) 
are composed of only two structural units corresponding to the two ^5 symmetric GBs. Molecular 
dynamics simulations are then conducted to investigate the mechanical response and the underlying 
deformation mechanisms of bicrystal models with different ^5 GBs under tension. Tensile deformation is 
applied under both 'free' and 'constrained' boundary conditions. Simulation results revealed different 
mechanical properties of the symmetric and asymmetric GBs and indicated that stress state can play an 
important role in the deformation mechanisms of nanocrystalline materials. 

It has long been recognized that grain boundaries (GBs) are important microstructure features and can 
significantly affect the properties of polycrystalline materials'. When the grain size is reduced to ultrafine 
or nano scale, the effects of GBs on the material properties become more significant since the traditional 
deformation mechanisms based on nucleation and propogation of lattice dislocation are replaced gradually by GB 
mediated processes, such as GB sliding"", grain rotation"* *, diffusional creep^, dislocation nucleation or absorption 
at GB etc.''"'. Much experimental work and many atomistic simulations have been conducted to examine various 
GBs on their energy, structure and properties. However, previous research work has focused primarily on 
symmetric GBs which possess mirror symmetry of crystallographic planes. In contrast, very few atomistic 
simulations have been conducted on the structure and related properties of asymmetric GBs, even though 
experimental observation has shown that most GBs in real polycrystalline materials are actually asymmetric"" 
and they can affect the material properties more significantly''''''. 

Recently, both experimental observations"" and atomistic simulations^^"^' have been attempted of asymmet- 
ric tilt GBs, and these give us a better understanding of the structures and energy on these boundaries and provide 
insight into the related GB properties, e.g., the structural transformations, dissociations, faceting transitions etc. 
Much of these works concentrate on the ^3 family because the most common ^3 GB is the coherent twin 
boundary with very low boundary energy and they are more frequently observed in polycrystals'""'"'^""'. Tschopp 
and McDowelP'' took advantage of atomistic simulation on both symmetric and asymmetric ^3 < 1 1 0 > tilt 
GBs and found that the structure and energy of the asymmetric GBs were closely related to the corresponding 
symmetric ones. Their investigation^' on dislocation nucleation from different ^3 asymmetric boundaries under 
uniaxial tension demonstrated that the properties of GB depend not only upon the misorientation between grains, 
but also upon the inclination of the GB plane. Lin^"" and Luo^"* performed simulations on the asymmetric ^3 < 1 1 
0>{1 10)/(1 14) GB in copper bicrystals under shock loading to examine the influence of GB on the mechanical 
behavior of the bicrystal system. Compared with the productive research findings on ^3 asymmetric boundaries, 
very few experiments or simulations have focused on other low index coincident site lattice (CSL) systems. 
Recently, Tschopp and coworkers'" investigated the structure and energy of several low CSL bicrystal systems, 
i.e. ^5, "^13 GBs around [0 0 1] misorientation axis and ^9, JZ^l around [1 1 0] misorientation axis. 
However, they did not study the correlation between the asymmetric GB structures and their mechanical 
behavior. H. Zhang performed a series of simulations to examine the influence of GB inclination on the mobility^\ 
diffusivity'''' and migration^*''^ of ^5 tilt asymmetric GBs, but no further research was conducted to analyze the 
mechanical response in relation to the deformation mechanisms at atomic scale. 



SCIENTIFIC REPORTS | 4:5919 | DOI: 1 0. 1 038/srep0591 9 



1 



In this study we carried out molecular dynamics (MD) simulations 
on two symmetric and four asymmetric ^5 GBs of copper under 
tensile deformation to probe their mechanical behavior and the 
related deformation mechanisms at atomic scale. Although the prim- 
ary aim of this work is to show the effect of GB inclination on tensile 
response, the influence of transverse stress on the mechanical beha- 
vior is also examined. In the simulations, the tensile deformation is 
conducted under both 'free' and 'constrained' tension boundary con- 
dition, which will be specified in the methods section. 

Results 

GB structure and energy. GB structures are identified using the 
common neighbor analysis (GNA) technique suggested by Schiotz 
et al.''''. Fig. 1 shows the equilibrium structures of the two symmetric 
and four asymmetric ^5 GBs of Cu at OK which were investigated. 
The black and white balls correspond to the two adjacent atom layers 
along the [0 0 1] tilt axis. Structure units as defined by Rittner" are 
used to illustrate the boundary structures. The six-member kite- 
shaped unit is referred to as the "E" unit which organizes the 
symmetric ^5 (3 1 0) ((f) = 0°) boundary as outlined in Fig. 1 (a). 
The symmetric ^5 (2 1 0) ((|) = 45°) is composed of the topologically 
identical structure unit "E"', which differs only in direction and 
arrangement along the boundary plane, as shown in Fig. 1 (f). 
Careful investigation of the four asymmetric GBs shows that all of 
them are organized by the combination of E and E' units that 
correspond to their symmetric ^5 (3 1 0) and ^5 (2 1 0) 
boundaries but with different ratios. In other words, the 
asymmetric GBs can be decomposed into their corresponding 
symmetric boundaries. For example, the period vector for the 
asymmetric (9 2 0)/(6 7 0) (4) = 30.96°) is ao [7 -6 0] or ao [2 



— 9 0] and this can be separated into one period of the ^5 (3 1 0) ((]) 
= 0°) symmetric GB and three period of the (2 1 0) ((j) = 45°) 
symmetric GBs. This reaction can be described by: 

ao[7-6 0]=ao[l -3 0]-h3ao[2-10] (l) 
orao[2 -9 0]=ao[-l -3 0]-|-3ao[l -2 0] (2) 

It should be noted that the ratio of E unit to E' unit decreases with 
the increase of inclination angle (j) and the equal ratio is found in the 
(1 0 0)/(4 3 0) ((j) = 18.43°) boundary, as shown in Fig. 1 (c).This 
finding that asymmetric GBs can facet into their corresponding sym- 
metric ones is similar to the finding of the 5^3 [1 1 0] grain boundary 
family^'. 

Since the structures of asymmetric GBs are closely related to their 
corresponding symmetric boundary structures, a faceting model is 
proposed to predict the energy of asymmetric GBs with idea faceting 
structures, i.e. the energy of the faceted asymmetric boundaries can 
be predicted by simply using a weighted fraction of the relative con- 
tribution of facet lengths times the respective symmetric boundary 
energy". This relationship can be expressed by the following equa- 
tion, 

r ^ .^/coscxx] [sin CD] , , 

Y, = Y,. ■ [cos CD - sm <D (— j J + Ys. [—J (3) 

where Ya is the predicted energy of the asymmetric GB with different 
inclination angles (j), Ysi and Ys2 are the calculated energy of the two 
corresponding symmetric GBs, a is the interval angle separating the 
two symmetric GBs which depends on the crystal symmetry around 
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Figure 1 | ^5 GB structures of Cu for various inclination angles (() at OK. The structures are viewed along the [0 0 1] tilt axis, atoms on consecutive 
(0 0 2) plane are shown as black and white. The boundary normal vector of grain A and grain B are marked on the right-hand side for each GB. The 
structure unit at each boundary plane are outlined by the solid line as marked by E and E'. (a) and (f) are the two symmetric GBs, (b)-(e) are the 
four asymmetric GBs. 
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Figure 2 | Grainboundary energy as a function of inclination angle of six J] 5 GBs. The solid square is the predict GB energy from equation (3), the dotted 
line is fitted by the predicted results which represents the trend of energy for ^5 asymmetric GBs. The solid dot is the GB energy of MD simulation. 



the tilt axis (a = 30" for [1 1 1] axis, cx = 45^ for [0 0 1] axis and a = 
90° for [1 1 0] axis). Here, the calculated energy of the two ^5 
symmetric GBs is 948 mj/m^ and 997 mj/m^ respectively, the pre- 
dicted energy value according to equation (3) and the calculated 
energy value from MD simulation of the four asymmetric GBs are 
plotted in Fig. 2. The consistent trend and small deviation between 
the predicted and the calculated energy values indicates that the 
faceting of asymmetric GBs into the structure units of the corres- 
ponding symmetric GBs is favorable from the point of view of energy. 

Mechanical response and deformation mechanism. Tensile 
simulations are conducted under free tension boundary condition 
and constrained tension boundary condition respectively, the 
simulation is specified in the method section. Fig. 3 (a) and (b) 
show the tensile stress-strain response of bicrystal models with 



different ^5 GBs at 10 K. For the mechanical properties, the 
system stress is attained by calculating the pressure of the entire 
system of atoms, system strain is derived from the positions of the 
periodic boundaries. On the whole, the maximum tensile stress of the 
two symmetric GBs ((|) = 0° and c|) = 45°) is higher than the peak 
stress of the four asymmetric GBs. The mechanical behavior of 
different ^5 GBs can be associated with their energy, as seen in 
Fig. 2. As calculated in the simulation, ^5 ((|) = 0° and (() = 45°) 
GB shows the comparative lower energy, so it may have the more 
stable structure and the tensile strength is higher. It is worth noting 
that the maximum tensile stress of each case under constrained 
boundary condition is significantly higher than the value of free 
boundary condition. The higher tensile stress can be attributed to 
the stress that developed transverse to the loading direction during 
the deformation process. 




Tensile Strain (%) Tensile Strain (%) 

Figure 3 | Mechanical response of 5]]5 GBs with different inclination angles (() at 10 K under (a) free tension boundary condition and (b) constrained 
tension boundary condition. 
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(a) e=0% (b) £=8.9% (c) e=9.04% 

Figure 4 | Snapshots of Cu bicrystal with J]]5 ((() = 0°) GB at different deformation stage under free tension boundary condition. Images are 
colored according to the CNA parameter. Atoms with perfect fee structures are removed to facilitate viewing of the defective structures. Atoms colored 
with yellow organize the GB plane and the dislocation core, while the blue atoms represent the stacking fault. 



Visual inspection of the MD simulation results indicates that the 
maximum tensile stress corresponds to the nucleation of partial dis- 
locations, in agreement with results for uniaxial tensile of Cu bicrys- 
tal by Spearot et al.-''". Fig. 4 shows the snapshot of atoms in Cu 
bicrystal with ^5 (4) = 0°) GB at different deformation stages under 
free tension boundary condition. Images are colored according to the 
CNA parameter. Atoms with perfect fee structures are removed to 
facilitate viewing of the defective structures. Atoms colored with 
yellow organize the GB plane and the dislocation core, while the blue 
atoms represent the stacking fault. The dislocation extraction algo- 
rithm (DXA)'" '"' is used to compute their Burgers vectors. The GB 
region becomes coarsening when the tensUe deformation is increas- 
ing until the maximum tensile stress has been reached. In Fig. 4 (b), at 
the beginning of the stress drop (e = 8.9%), the image shows that 
partial dislocation loops with both edge and screw character are 
nucleated from the bicrystal interface into the upper grain and the 
lower grain simultaneously. DXA analysis indicates that Shockley 
partial dislocation with Burger's vectors b = (1/6)[1 1 2] and b = 
(1/6)[1 1 —2] nucleated from the bicrystal interface and slip on the 
(111) and (1 1 — 1) plane. According to the Schmid factor analysis, 
they are the most favorable slip system. The tensile stress required to 
nucleate the first partial dislocation from the 5^5 (c() = 0°) GB at 10 K 
is calculated as 6.28 GPa, which corresponds to a critical resolved 



shear stress of approximately 3.08 GPa for the given lattice orienta- 
tion. The result is comparable to the maximum resolved shear stress 
under the uniaxial tensile deformation of bicrystal Cu at 300 K by 
Spearot"". With the increase in the tensile strain, dislocations nucle- 
ate continuously from the GB plane and slip in each grain, as seen in 
Fig. 4 (c). 

Fig. 5 shows the atomic details of Cu bicrystal with ^5 ((f) = 
11. 31°) GB at different deformation stages under free boundary con- 
dition. Atoms are colored by the CNA parameter and remove the 
atoms with fee structures to facilitate the defective structures. The GB 
region expands and becomes coarsened when deformation occurs. 
Then, partial dislocations with Burger's vectors b = (1/6) [1 1 2] and 
b = (l/6)[ll— 2] are nucleated in the lower grain region, as shown in 
Fig. 5 (b) at 8 = 7.3%. However, unlike the case of (4) = 0°) GB, 
dislocation nucleates into only one crystal lattice when the maximum 
tensile stress has been reached. This phenomenon can be attributed 
to the asymmetric GB with different orientation angles in the two 
lattices. Partial dislocations are nucleated and emitted continuously 
into the lower grain until the tensile strain reaches s = 7.6%, as 
shown in Fig. 5 (c). The slip system is now activated in the upper 
grain, evidenced by a partial dislocation with Burger's vectors b = ( 1/ 
6)[1 1 2] nucleated from the interface and slip along the (1 1 — 1) 
plane. After that, dislocation slips collectively in both grain regions 





(a)E=0% 



(b) £=7.3% 





(c) £=7.6% 



(d) £=7.78% 



Figure 5 | Snapshots of Cu bicrystal with 5]]5 ((() = 11. 31°) GB at different deformation stage under free tension boundary condition. Images are 
colored according to the CNA parameter. Atoms with perfect fee structures are removed to facilitate viewing of the defective structures. Atoms colored 
with yellow organize the GB plane and the dislocation core, while the blue atoms represent the stacking fault. 
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accommodate the plastic deformation during the tensile process. The 
same phenomena has been observed in other cases of Cu bicrystals 
with asymmetric GB, as shown in S-Fig. 2 

Fig. 6 shows the different deformation stages of Cu bicrystal with 
^5 ((|) = 18.43°) GB when subjecting to tensile deformation under 
constrained boundary condition. Fig. 6 (a)-(c) gives the projected 
view of the atomic snapshots, atoms with perfect fee structure are 
rendered as dark blue, atoms at GB area and the dislocation core are 
rendered as yellow, the light blue atoms indicate the stacking fault. 
Dislocations are extracted and they are converted to continuous 
lines, as shown in Fig. 6 (d)-(f). In the elastic deformation stage, 
all the boundary atoms are well organized and the GB keeps its 
equilibrium configuration until the maximum tensile stress is 
achieved. Fig. 6 (a) corresponds to the initial stage of the plastic 
deformation when a crack is initiated on 'grain boundary 1' and 
'grain boundary 2' simultaneously (as introduced in Methods section, 
the periodic boundary condition applied in the Y direction intro- 
duces a second boundary plane). Atoms are beginning to shuffle at 
the crack tips where partial dislocation has nucleated, as seen in the 
enlarged area in Fig. 6 (a). Again, DXA analysis indicates the slip of 
partial dislocations occurring on both activated (111) and (1 1 —1) 
planes in the lower grain region, which is in agreement with Schmid 
factor analysis. As the tensile strain increases, the fast decline of 
tensile stress is associated with the progressive separation of the 
GBs as well as the partial dislocation nucleation from the crack tips, 
as seen in Fig. 6 (b) and (e). The high density of the dislocation 



network at 8 = 8.44% only stays for a short period and soon decreases 
gradually with the cleavage of the GB plane and the separation of the 
two grains, as seen in Fig. 6 (c) and (f). Interestingly, the dislocation is 
not seen to be emitted into the upper grain region during the tension 
process. Unlike the massive dislocations and their slipping in the 
bicrystal models under free tension boundary condition, only limited 
dislocations are observed in the cases studied under constrained 
boundary condition. The crack initiates and extends along the GB 
plane in a cleavage manner and eventually results in the fracture of 
the model, as shown in S-Fig. 3. 

Discussion 

Different from the symmetrical GBs which possess mirror symmetry 
of crystallographic planes, asymmetric GBs are relative complex with 
a multiplicity of atomic structures resulting from the boundary dis- 
sociation, nano-faceting, and other fundamentally interesting struc- 
tural effects. Simulations^''''" revealed that asymmetric GBs tend to 
break into nano-scale facets which are composed of their corres- 
ponding symmetric boundaries. Sutton and Balluffi*" proved that 
the faceting of an asymmetric GB into a symmetric GB is a geomet- 
rical possibility. This is also confirmed in our simulation of the ^5 
GBs where all the asymmetric boundaries are organized by the com- 
bination but different ratio of E and E' structure units which corre- 
spond to the 4) = 0° and (|) = 45° symmetric boundary respectively. 
However, a recent study of JI^H [1 1 0] tUt GBs*' with different 
inclination angles revealed that it is not necessary for the asymmetric 
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(d) 8=8.34% 



(e) 8=8.44% 



(f) 8=8.8% 



Figure 6 | Snapshots of Cu bicrystal with "Y^S ((() = 18.43°) GB at different deformation stage under constrained tension boundary condition. 

Images (a)-(c) are colored according to the CNA parameter, atoms with perfect fee structure are rendered as dark blue, atoms at GB area and the 
dislocation core are rendered as yellow, the light blue atoms indicate the stacking fault. Dislocation segments are extracted in (d)-(f). 
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boundaries to dissociate into symmetric ones. The faceted boundary 
may do not even belong to any particular CSL. It was reported that 
the asymmetric "^11 [1 1 0] GBs dissociated into a low- angle GB 
formed by Shockley partial dislocations and a high-angle non-^11 
boundary. This kind of faceting into non-symmetric boundary struc- 
tures and incorporating facets of non-^ boundary are not found in 
our simulation of ^5 GBs, primarily because of the differences in 
boundary energy. Gokon et al.''''''''' experimentally determined the 
boundary energies of [1 1 0] and [1 1 0] tilt GBs in Cu 
for various inclination angles and, found the facet planes often corre- 
spond to a nearby symmetric or asymmetric boundary with low 
energy. According to the simulation result of ^5 GB energy in 
our work, as seen in Fig. 2, the energy of asymmetric GBs is com- 
parable with the 4) = 45° symmetric boundary and higher than the c|) 
= 0° symmetric boundary. From the point of view of energy, this 
may explain why such faceting of ^5 asymmetric GBs occurs intend 
to the symmetric (|) = 0° and ^ = 45° planes. 

Simulations reveal that dislocations are nucleated and emitted into 
both grains of the symmetric GBs once the maximum tensile stress 
has been reached, while this emission occurs only at the lower grain 
region in the asymmetric GBs at the beginning, and then the slip 
systems in the upper grain can be activated, as shown in S-Fig. 2. This 
phenomenon is expected for asymmetric GBs with inclination 
angles. Specifically, grains on each side of an asymmetric GB are 
oriented differently relative to the tensile stress direction, resulting 
in different Schmid factors. Therefore, the slip systems can be acti- 
vated easily in the grain which associates with the higher Schmid 
factor, since the slip systems in this grain have higher resolved shear 
stress. Fig. 7 shows the maximum Schmid factor for both grains of the 
bicrystal models as a fiinction of the inclination angle, from which we 
can see that the maximum Schmid factor is higher in the lower grain 
in all of the four cases of asymmetric GBs. This may explain why the 
partial dislocations preferentially nucleate first into the lower grain. 
For the asymmetric GBs, the nucleation of dislocations in the upper 
grain only occurs at higher strains when excessive dislocation emis- 
sion severely alters the initial orientation of the lower grain. 



In this study, we found that stress state can play an important role 
in the deformation mechanisms of nanocrystalline materials. 
Specifically, massive dislocation nucleation and propagation is the 
dominant mechanism during the tensile deformation under free 
tension boundary condition, as shown in S-Fig. 2. However, due to 
the transverse stress that applied perpendicular to the tensile stress 
under constrained tension boundary condition, the nucleation of 
dislocation and its propagation has been restrained. Meanwhile 
cracks are prone to initiate and can extend rapidly along the GB 
plane, which becomes the dominant mechanism during tension 
and leads to a more brittle mode of failure, as shown in S-Fig. 3. 
Moreover, only partial dislocation loops are nucleated from the Cu 
bicrystal interfaces during the tensile process under free tension 
boundary condition. However, it is interesting to find that full dis- 
location loops can also nucleate from boundary plane under the 
constrained tension boundary condition, as shown in Fig. 8 for the 
case of ^5 ((|) = 0°) GB. I Fig. 8(a), the leading partial dislocations 
have nucleated and moved away from the boundary plane, leaving 
behind an intrinsic stacking fault. As deformation proceeds, the 
trailing partial dislocation is evidenced to emit, resulting in a fuU 
dislocation loop and then pass through the periodic boundaries, as 
shown in Fig. 8(b) and (c). In previous reports*'* the trailing partial 
dislocation is considered can mainly nucleate from the fee metals 
with high stacking fault energy {e.g. Al 146 mj/m^"*^), while the nuc- 
leation of fuU dislocation in the fee metals with low stacking fault 
energy {e.g. Cu 44.4 mj/m'*'') is thought to be difficult. To our under- 
standing, few of the previous simulation studies has reported the full 
dislocations nucleation from Cu sample. Result of this study that full 
dislocations can be nucleated and propagate in Cu bicrystals under 
constrained tension boundary condition in return emphasizes the 
important role of stress state in the deformation mechanisms of 
nanocrystalline materials. 

Methods 

Model construction. In this study, bicrystal model with fixed orientation of the GB 
plane (0 — 36.9' for ^5 GBs) and fixed tilt axis ([0 0 1] axis) is employed, which 
enables a more controlled investigation of specific GB properties. Fig. 9 shows a 
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Figure 8 | Snapshots of Cu bicrystal with ^5 = 0°) GB at different deformation stage under constrained tension boundary condition. Images 
are colored according to the CNA parameter. Atoms with perfect fee structures are removed to faciUtate viewing of the defective structures. Atoms 
colored with yellow organize the GB plane and the dislocation core, while the blue atoms represent the stacking fault. 



schematic of the computational cell used in our simulations. A bicrystal model is 
created by constructing two separate crystal lattices (grain A and grain B in Fig. 9) 
with different crystallographic orientation and joining them together along the Y axis. 
The misorientation is defined as angle 0 between the [10 0] direction of the two single 
crystal grains, and the inclination is defined as angle (f) between the bisector of the 
misorientation and the boundary plane. The values of (f) considered in this study and 
the Miller indices of the boundary plane are presented in S-Table 1. Due to the 
fourfold symmetry of the fee lattice, the inclination angles from 0° to 45° cover all 
distinct boundary structures. 



Grain A 



Grain B 




Figure 9 | Schematic of Cu bicrystal with a ^5 [0 0 1] asymmetric tilt 
grain boundary. The misorientation is defined as angle 0 between the [ 1 0 
0] direction of the two single crystal grains and the inclination is defined as 
angle <\> between the boundary plane and the bisector of misorientation 9. 



Periodic boundary conditions are used in all directions (X, Y and Z). Specifically, 
periodic boundary condition are applied to X and Z directions to simulate an infinite 
boundary plane between the two grains which can eliminate the effect of free surface. 
The periodic boundary condition in Y direction introduces a second boundary plane 
into the model. For each initial configuration, it is important to adjust the model size 
to construct identical atomic structures of the two boundary planes to ensure they 
have the same equilibrium structures and energies after energy minimization. 
Otherwise, the different equilibrium structures of the two boundary planes will result 
in, for example, a higher energy metastable plane at grain boundary 1 and a stable 
plane at grain boundary 2, which can significantly affect further simulation results. A 
number of initial "starting positions" of grain A and grain B are tested to find the 
possible GB structures'"^ Molecular statics calculations which employ a nonlinear 
conjugate gradient algorithm are conducted on all the tested GB structures to 
determine their minimum energy configurations. 

GB energy calculation. After the procedure of energy minimization, the energy of 
each equilibrium GB structure is then calculated and compared to find the possible 
global minimum energy configuration. In this study, The GB structure with the lowest 
energy is regarded as the stable structure, while other GB structures with higher GB 
energy are reckoned to be metastable which will not be considered in our study. Since 
the periodic boundary condition in Y direction generates a pair of boundary planes in 
the bicrystal model, the energy associated with the grain boundary is calculated by the 
following equation: 



Tgb - 



2A 



(4) 



where Esystem is the system energy of the equilibrium bicrystal model, Eatom is the 
potential energy of a single atom in the perfect Cu lattice ( — 3.54 eV), N is the total 
number of atoms contained in the model and A is the area of the GB plane. 

Molecular dynamics simulation. All the simulations in this study are performed 
with the parallel molecular dynamics code LAMMPS'^^ The embedded-atom method 
(EAM) potential developed by Mishin et al. for copper^^ is employed, which can fit a 
large set of experimental and first- principles data. After the minimum energy 
configuration is attained, the simulation model is equilibrated using MD in the 
isobaric-isothermal (NPT) ensemble at a pressure of 0 bar and a temperature of 10 K 
for 20 ps. A constant rate of 2 X 10^ s~^ is applied perpendicular to the boundary 
plane (along Y direction) at a temperature of 10 K. Tensile deformation is performed 
under either 'free' or 'constrained' boundary condition. These boundary conditions 
are very similar to those used by Kitamura et al." and Spearot et al.'^' to study the effect 
of boundary condition on tensile deformation of nickel single crystal and copper 
bicrystal respectively. Under free tension boundary condition, the boundaries in the 
lateral directions are allowed to expand or contract during the deformation process 
and the transverse stresses are kept free (Qxx ~ ^zz — 0)- Under constraint tension 
boundary condition, computational models are strained at a constant rate along the Y 
axis while keeping the model dimensions along the X and Z axis fixed (e^x ~ £zz — 0)- 
This boundary condition considers the transverse stress along the X and Z axis during 
the tensile deformation process. The schematic of tension simulation is shown in S- 
Fig. 1. The visualization tools Atomeye^^ and Ovito^' are used to illustrate of the 
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bicrystal models. The common neighbor analysis (CNA) technique'^ is used to 
identify the defect structure and its evolution during the simulations. The dislocation 
extraction algorithm (DXA)^'^^^ is used to convert identified dislocations into 
continuous lines and compute their Burgers vectors. 
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